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O ' ABSTRACT 
<N 

^ '. We propose a modified version of the X-ray spectral index and an intrinsic 

cutoff frequency of inverse Compton radiation from the brightest knot of the 



> 



c3 



CO '. M87 jet, in conjunction with an application of the new conceptions of injection 

and diffusive shock acceleration (DSA) of electrons in magnetized filamentary 
plasma to the specified source. The drop of the X-ray flux density in a transitive 
frequency region is associated with the interplay of ordinary synchrotron cooling 
and weaker magnetic fields concomitant with the smaller scale filaments that 
allow the electron injection, while the radio-optical synchrotron continuum is 
dominantly established by the major electrons that are quasi-secularly bound to 
larger filaments. With reference to, particularly, the updated external Compton 
model, we demonstrate that in the Klein-Nishina regime fading inverse Comp- 
tonization, the injected electrons can be stochastically energized up to a Lorentz 
factor as high as 5 x 10 10 in the temporal competition with diffuse synchrotron 
cooling; this value is larger than that attainable for a simple DSA scenario based 
on the resonant scattering diffusion of the gyrating electrons bound to a sup- 
posed magnetic field homogeneously pervading the entire knot. The upper limits 
of the photon frequency boosted via conceivable inverse Compton processes are 
predicted to be of the common order of ~ 10 30 Hz. The variability of the broad- 
band spectrum is also discussed in comparison to the features of a blazar light 
curve. The present scenario of a peta-eV (PeV; 10 15 eV) electron accelerator, the 
"Pevatron," might provide some guidance for exploring untrod hard X-ray and 
gamma-ray bands in forthcoming observations. 

Subject headings: acceleration of particles — galaxies: individual (M87) — galax- 
ies: jets — gamma rays: theory — magnetic fields — radiation mechanisms: 
nonthermal 
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INTRODUCTION 



A fundamental question — the observable highest frequency of radiation from extra- 
galactic jets — remains unresolved, as the closely relevant physics of the particle acceleration 
incorporated with injection and the cosmological transport of energetic photons is largely 
unknown in both theoretical and observational aspects. As a candidate celestial object to 
study those issues in detail, of a suitable nature is a nearby powerful radio source, par- 
ticularly the M87 jet, and so far, many eff orts have been devoted to measure the photon 
spectr a in higher frequency bands (optical, IMeisenheimer et al.lll996l ; iPerlman et al.l 12001 
X-ray 



i m mgner irequency Danas ( optical, livieisenneimer et ai-lli^Dl ; irermian et ai-llzuuii ; 
Marshall et al-l booj Iwilson fc Yanell2002h. In the previous key work that was based 



on ground-based observations (IMeisenheimer et al.lll996l ). the synchrotron spectra were sug- 
gest ed to, in the entire jet, have an almost consta nt cutoff of (0.8 — 2.5) x 10 15 Hz (see 



also Perez- 



burnon et al.l Il988l ; IStiavelli et al.l 119971 ) , compared to a universal value among 



radio jets ( Keell Il988l ). The corresponding Lorentz factor of the electrons was estimated 
to be 9 x 10 5 , which could be deduced from a simple model of the diffusive shock accel- 
eration (T2SAJ_th^ 1 t_j^hed_oji_t^ scattering by the turbulent Alfven 
waves (IBiermann fc Strittmatterl 119871 ; iLongairl Il992l ). However, state-of-the-art Chandra 
observations revealed that the spectra indicate no steep cutoffs, extending to X-ray regions 
( Marshall et al. I l2002h . far beyond the previous cutoff. The X-ray continua have been ascribed 
to synchrotron radiation, whereas inverse Compton (IC) models as a dominant X-ray mech- 
anism have been ruled out. The Lorentz factor of the X-ray-emitting electrons is considered 
to be of order ~ 10 6 — 10 8 (IMeisenheimer et al.lll996l ; IPerlman et al.ll200ll ; iMarshall et al. 



20021 ; IWilson fc Yangj|2002l ). although the achievable maximum is as yet poorly known. 



A clue to this puzzling question is that simple models invoking shocks an d synchrotron 
losses in a uniformly filled j et overpredict the X-ray flux density by large factors feiretta et al. 



199ll ; IMarshall et al.ll2002l ). In fact, the fine structure of magnetized filame nts including; the 



limb-brightening feature has been found at many positions along the jet (IReid et al.l 11989 



Owen et al.l Il989l ; IStiavelli et al.l Il99ll ; lLobanov et al.l 120031 ) . The cancellation of the small- 
scale magnetic fields could be responsible for the observational fact that in the interior of 
the jet, the magneti c field is significantly free (lOwen et al.lll989l ). reducing the polarization 
( ICapetti et al.lll997l ). Accordingly, the requirement for modeling a nonuniform jet is that in 
the interior the turbulent state ought t o be rather vital, in contrast to conventional jet models 
involving a turbulent radial edge (e.g. 



Ostrowskill2000l ) caused by surface-mode instabilities 



such as the Kelvin-Helmholtz instability. 



At this juncture, the novel model by lHonda fc Hondal (120021 ). in which the jet is hypothe- 
sized to co mprise the magnetize d filaments that are self-organized via a current filamentation 
instability (IHonda et al.ll2000al ). seems to satisfy this requirement. In their model, the vec- 
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tors of fluctuating magnetic fields in the interior are supposed to be randomly oriented in 
the transverse direction of the filaments, such that they could be macroscopically canceled 
out, although in the envelope r egion, the uncanceled, ordered magn etic field could ap pear 



( IHonda fc Hondall2004al . l2005al ). By invoking a similar environment, iFleishmanl (120061 ) has 
recently examined synchrotron radiation from the electrons injected into the small-scale 
random magnetic fields and found that the diffuse emission including higher harmonics re- 
produced the observational X-ray continuum at the knots B and CI. However, the mecha- 
nisms of the electron inj ection reflecting the ene rgy hierarchy that ought to be involved in 



the filamentary plasma (IHonda fc Honda! l2005al ) and the related in situ acceleration were 



not self-consistently taken into consideration. It is important to investigate these details, 
because the complex interplay will actually connect the confirmed radio-optical continuum 
with the diffuse synchrotron spectrum in th e energetic regime to give a transitive spectrum 



compared to the observed X-ray spectrum ([Marshall et al.ll2002l ). which is expected to de- 
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their trajectory in the filamentary (stochastic) medium and thereby have a shorter mean free 
path (mfp) that deviates from the universal gyroscaling (valid for the simple case in which 
a mean magnetic field can be well defined); this is just the case that is advantageous to the 
DSA. The possible reboost of electron energy will, in turn, upshift the IC gamma-ray cutoff 
frequency. 

In the context of the filamentary jet model, here we propose a novel scheme to con- 
struct the steepened X-ray synchrotron spectrum (§ 2) and attempt to provide genuine 
uppermost energies of an accelerated electron and photon at knot A (§ 3), where X-ray flux 



is so prominent as to be comparable to that from the galactic core (IMarshall et al.l 12002 



Perlman fc Wilsonl 120051 ). The theory is based on the DSA model incorporated with the 
new type of injection reflecting electron escape from small-scale magnetic traps, which is 
referred to as the "transition injection" (§ 2.1). We obtain a modified X-ray spectral in- 
dex consistent with the large flux density reduction measured in the Chandra regime and a 
nonlinearly evolving feature of the current filamentation insta bility (§ 2.2). In the conceiv- 



able radiation fields, including external Compton (EC) targets (jStawarz et al.ll2003l ). we find 
that the injected electrons can be, in situ, accelerated up to the 10 PeV range (§ 3.1). It 
is shown that the flux density of IC radiation surpasses, above the soft gamma-ray region, 
that of synchrotron radiation, and the highest frequency reaches the theoretical upper limit 
reflecting the maximum electron energy (§ 3.2). The IC spectrum can be accommodated 
with the recent High Energy Stereoscopi c System (HESS) o bservation of sub-TeV gamma 
rays emanating from the galaxy region (IBeilicke et al.l 120051 1. although the origin of these 
gamma rays is still unresolved. The energy budget of the diffuse synchrotron radiation in 
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the gamma-ray spectral tail is also estimated (§ 3.3). Moreover, we include a discussion con- 
cerning the variability of the broadband spectrum (§ 4). The presence of fine filaments turns 
out to be essential for the transition injection, raising the flux density level of the diffuse 
radiation and engendering the high-frequency variation, and therefore the observability of 
such filaments is also discussed (§ 5). 



2. SYNCHROTRON SPECTRUM EXTENDING TO THE X-RAY BAND 

WITHOUT STEEP CUTOFF 

Highlighting the conception of the transition injection that mediates the energy hierar- 
chy involved in a magnetized filamentary plasma, we consider the electron preacceleration 
below the injection energy. First we recover the radio-optical synchrotron continuum in the 
cold limit (§ 2.1), and second, we construct the transitive X-ray spectrum (§ 2.2). 



2.1. Impacts of the Small Spatioscale Entailed in the Transition Injection 



Begin by considering a low-energy state in which cosmic-ray (test) electrons are bound 
to the magnetized current filaments with various transverse sizes A, exhibiting gyromotion. 
Apparently, this can be compared to the low-/3 plasma state, and along with this, it seems to 
be reasonable to suppose that the locally ordered fields, with a strength of |B(r)|, permeating 
through the filaments, are preferentially disturbed by Alfven waves. As is accommodated 
by the test particle approximation, we presume a small-amplitude disturbance that does 
not affect the zeroth-order magnetic fields sustained by the bulk filaments (the validity is 
confirmed later). In this aspect, recalling the consensus that the jet knots are associated 
with shocks, the gyrating cold electrons can presumably be accelerated by a conventional 
Fermi-based DSA mechanism relying on the weak Alfvenic fluctuations. For a possible case 
of the scattering; center compre ssion ratio of r ~ 4 (reflecting a high Mach number shock; 
Vainiol Il999l : ISchlickeiseri booj ). the predicted synchrotron emission from the nonthermal 
electrons is certainly compatible with the observed radio-optical continu um of the power-law 



flux density distribution [F v oc v~ a \ which has indices around a ~ 0.5 (IBiretta et al.l 1991); 
that is, assigning the a-value to the canonical expression ct\ = (p — l)/2 (e.g., Longaii[l994 ) 
gives a differential spectral index of the electron density of p = (r + 2)/(r — 1) ~ 2 (for 
dn/de = ne~ p , where e = jmc 2 ). It is precisely noted that in general, the relation of r to 
the gas compres s ion ratio depends on the plasma conditions, such as the plasma /5-value 



Vainiol Il999l ; ISchlickeiseri 120021 . for the detailed calculations). As long as the gyroradius 



(see 

of the accelerated electrons r g>e is smaller than the transverse size of the local filament 
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trapping them (~ A), such a DSA mechanism will safely operate. With regard to energy 
restrictions, the acceleration of the electrons bound to a large filament ought to be knocked 
down by radiative losses before the escape from the filament (temporal limit), whereas the 
electrons bound to a small filament could escape via the energization (spatial limit). It 
is, thus, anticipated that there exists a critical scale of filament, A c , that discriminates 
between the large and small scales. Furthermore, with reference to the consequence of 
kinetic simulations that the amplification of magnetic fields a ccompanies the coalescence of 
filaments accumulating larger filaments JHonda et aliboOOal lbl). we consider a situation in 
which the filaments with larger (smaller) values of A involve the larger (smaller) values of 
|B|. Here we invoke the scaling of |B| oc A^ -1 ^ 2 , where (3 corresponds to the power-law 
spectral index of the magnetized fila mentary turbulence (possibly larger than the spec tral 
index for the Kolmogorov turbulence; iMontgomery fc Liulll979l ; iHonda fc Hondall2005al ). 



In order to quantitatively examine the A dependence of the maximum Lorentz factor of 
an accelerated electron 7*, of practical use is an extended scheme that can treat the oblique 
mean mag netic field with respect to the shock normal direction (for the methodological 
details, see lHonda &: Hondall2004bl ). We consider a simple model in which the energy density 
ratio of the fluctuating magnetic field to |B| 2 / (8tc), denoted as b (^C 1), is spatially constant, 
and the scattering of electrons occurs when the wave-particle resonance condition of r 9 e ~ 
k~ l is satisfied, where the wavenumber k is supposed to be in the inertial range of the 
power-law fluctua tion spectrum with an index o f 0' (e.g., 0' = 5/3 for the Kolmogorov- 



type fluctuations; iBiermann fc Strittmatterl 119871 ; IHonda fc Honda l2004bl ). Note that the 
local minimum of k (i.e., the correlation length of the Alfvenic fluctuations) is compared to 
the filament size A, whose maximum is described as ID, where / (< 1) and D denote the 
reduction factor of the filament correlation le ngth and the diame ter of the jet, respectively 
(typically, ID = 110 pc is expected at knot A; iBiretta et al.lll99ll ). 



As concerns the radiative losses of accelerated electrons, we take normal synchrotron 
and IC cooling into account. The former is ascribed to the gyration in the local mean field, 
whose intensity can be expressed in the phenomenological form of |B| = B m (\/ID)^~ 1 ^ 2 , 
where B m = |B| A=i - D . The diffuse effects considered in § 3 are ignored because b <C 1. For 
the latter, it should be remarked that while the synchrotron self-Compton (SSC) process can 
be preferentially ignored (§ 3.1), influential EC candidates do exist , including the starlight 



and extended dust emission of the host gal axy (IStawarz et a. 



a distance of 1 kpc from the galactic core (IHonda fc Honda 



2003). At knot A, which is at 



2004al ). the energy densities of 



the target photons are, in the comoving (shock rest) frame, suggested to be of the levels of 
M p h ~ 10~ 9 r 2 ergs cm -3 (stars) and ~ 10~ n r 2 ergs cm -3 (dust), where T = (1 — U 2 /c 2 ) -1 / 2 
is the Lorentz factor of the shock. The cosmic microwave background (CMB), which arranges 
u p h ~ 4 x 10~ 13 T 2 ergs cm -3 , plays a minor role in determining 7*. These EC scatterers 
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would uniformly pervade at the concerned region, such that u p ^ 
contrast to the inhomogeneous |B| (which is dependent on A). 



const is assumed, in 



In Figure 1, for a dominant stellar emission level of w p h = 10~ 9 ergs cm -3 and a promis- 
ing parameter set of /3 — 4.3 (for the reasoning, see § 2.2), j3' = 5/3, and ID = 110 pc, 
we show the A dependence of 7* for a fixed b = 10~ 3 , given B m = 100 fiG — 10 mG as a 
parameter (top), and for B m = 1 mG, given b = 10™ 4 — 10~ 2 (middle). Here we have set the 
compressio n ratio to r = 3.2 (f or p ~ 2.3, inferred from the radio-optical index of a = 0.67 

, the shock speed to U — 0.5c (as compared to the pattern 



at knot A; iPerlman et al. 



2001 



Biretta et al. 



19951). and the inclination angle of the local mean field to the 



speeds at knot A; 

maxi mum: cos -1 (U/c) = 60° (e.g.. lOwen et al.lll989l ; ICapetti et al.lll997l ; iHeinz fc Begelman 



19971 ). It is found that in any case, 7* appears to take the peak value at a certain value of 
A = A c , and as seen in Figure 1 (top), both 7* = 7*|a=a c and A c decrease as B m increases. For 
the electrons bound to the small filaments with a size of A < A c , the 7* value is determined 
by the simple relation of r gje (cx. 7*/|B|) ~ A that yields 7* oc A^ 1 ^ 2 , while for the electrons 



bound to the large filament with A > A c , a radiative loss limits the acceleration. For the 
former spatial limit, the 7*-value is independent of b, as seen in Figure 1 (middle). As for the 
latter temporal limit, in general, the balance betwe en the mean acceleration time, which has 
an approximate dependence of oc r„ 

,e(^/ r g,eY'~~ 1 flHonda fc Hondal l2004bh . and either the 
synchrotron cooling time, oc (7*|B| 2 ) -1 or the IC time, which has dependencies of oc (7*) _1 
(in the Thomson limit) and (hi7*) _1 (in the high-energy limit of the Klein-Nishina regime), 
yields the scaling of the form of 7* oc A -M , where 




= [(/? + 1) /?' - 2] / [2 (3-/3')], 
[(/3 + 1) (/?' - 1) - (/? - 1)] / [2 (3 
^ T [(3-/3')/ (2 -/3% 



respectively. In the displayed cases, the shortest loss time is that of the synchrotron cooling, 
such that the scaling with /i = /i syn = 2.6 appears, although for a less likely case in which 
both B m and b are small, an IC limit with [i ~ /xkn — 0.35 can appear (not shown in the 
figure). 

In Figure 1 (bottom), for the convincing param eters (at knot A) of B m = 4.59 mG (equal 



to the upper limit suggested by lOwen et al.lll989l ) and b = 10 3 (see § 2.2), we explicitly 



show the discrimination between the spatial (dotted line) and temporal (solid curve) limits, 
indicating the DSA tracks for gyrating electrons trapped in the filament with A < A c = 1.7 pc 
(solid arrows) and A > A c (dotted arrows). It is also demonstrated that in a high-energy 
regime and the possible frequency range of the target photon field (see § 3.1), the Klein- 
Nishina effects reduce the IC loss, fading the Thomson regime with /ix = 0.087 (thin dot- 
dashed line); as a result, the synchrotron cooling solely contributes to determine 7* for the 
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temporal limit regime, yielding a single scaling with /i syn (= 2.6, without /x-flip; soM curve). 
In any case, one can identify a critical Lorentz factor of 7* = 7.6 x 10 9 at which all minor 
electrons bound to the small filaments with A < A c can escape from them ( "escape tracks" ) 
to spaciously meander in a forest of the magnetized filaments. For the escaping electrons, 
more efficient energization is anticipated, as explained in § 3.1, wherein the critical energy 
7*mc 2 is pertinently compared to the injection energy for the additional DSA. We note that 
the major electrons quasi-secularly bound to the large filaments with A > A c ("secular-trap 
tracks" ) cannot reach the critical level of 7* because of the synchrotron cooling. 

This type of injection associated with the escape from the gy rotraps (namely, the bou nd- 



free transition) has been referred to as the transition injection (IHonda fc Honda l2005al ). It 
should be again noticed that the "escape" does not mean the leak of energetic electrons from 
the jet, but merely the escape from the magnetized fine filaments with A < A c ; that is, for 
the escape electrons to be additionally accelerated, of course, they still need to be confined 
in the shock accelerator. More details of the additional DSA are discussed in 5 3.1. 



2.2. X-Ray Spectral Index Linked to Turbulent Spectral Indices 



In the low-energy regime below 7*, the gyrating electrons bound to the filaments con- 
tinuously emit polarized synchrotron photons in an ordinary fashion. Thus, for the con- 
ventional Fermi-based DSA that engenders the power-law electron distribution (denoted in 
§ 2.1), the resulting flux density of syn chrotron radiation conforms to the generic scaling 



F v oc kv ( p 1 )/ 2 |B|^ +1 ^ //2 ( Longairlll994j ). For simplicity, we set k = const at the moment. If 



we recall the scaling of |B| oc A^" 1 ^ 2 , then it follows that the most dominant contribution 
to F v comes from the electrons bound to the filaments with an outer scale of A ~ ID. On 
the other hand, the familiar scaling of the synchrotron cutoff above which F v exponentially 
decays, v* oc (7*) 2 |B| u* oc 7* 2 |B| (valid for 7* > 1), is found to lead to v* oc A^ 2 ^^ 1 ^ in 
the secular-trap tracks of A c < A < ID (cf. Fig. 1, bottom). Hence, for values of \x > (j3 — 1)/4, 
which are typical for fi = yU syn (/5 / 7^ 1), the F v spectrum includes the most dominant com- 
ponent having a lowermost u*, which is considered to dominantly contribute to establish 
the so-called "break frequency" (commo nly identified in the optical to soft X-ray band; 
Perlman et al.ll200ll ; iMarshall et al.ll2002l ). The property that the less dominant component 
involves a higher u* cooperates to steepen the F v spectrum, in accordance with the X-ray 
measurement that indicates that the F v level falls below that inferred from extrapolation of 
the featureless radio-optical continuum to the X-ray region. The F v drop in the frequency 
change from v\, = v* 



\=ID 



to v* = v* 



A=A C 



can be estimated as 

(p+l)09-l)/4 



Xc/(ID) 



(2) 
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Note that equation (jSJ) is, through A c , dependent on B m and b (see Fig. 1, top and middle) 

In Figure 2, for (/?, (3') = (4.3, 5/3), given b = KT 4 - KT 2 , we plot u b (top), u* (middle), 
and 4> (bottom) as a function of £? m for the typical value of ID = 110 pc (thick curves) and 
a shorter value of ID = 11 pc (thin curves). The shock conditions have been chosen to 
be the same as those for Figure 1. It appears in Figure 2 (top) that u b decreases as B m 
increases, according to the dependence of u b oc i? m ' 4Alsyn ^ /3 l ' x ' = -B~ 21 in the regime in 
which the synchrotron cooling definitely limits the value of 7. In the major ranges of B m 
and b, as seen in the middle panel, u* remains almost constant, since the change of |B|a=a c is 
canceled out by that of (7*) 2 in the synchrotron-limiting regime. Also, in the corresponding 
regime, along the dependence of <fi oc A 2 ' 8 , the decrease of A c with increasing B m (Fig. 1, 
top) leads to the decrease of (j), as depicted in the bottom panel of Figure 2. Note that for 
a small value of B m (leading to a small | B | a=a c ) 5 the appearance of the synchrotron limit 
can be incomplete for a small b that lowers 7* (Fig. 1, middle) to weaken the Klein-Nishina 
effect; this property yields the peculiar B m dependencies of u* and in the marginal (B m , 
6)-domain, as shown in the middle and bottom panels of Figure 2. The values of u b , u*, and 
cf) are likely enhanced for large values of b, and, in addition, u b and <fi are further boosted 
for the small value of ID, although u* is unchanged, except for the IC-limiting regime. 
With regard to a comparison wit h the observational data, for a confirmed radio-optical F u 
spectrum with index a = 0.67 (IBiretta et al.l Il99ll ; iPerlman et al.ll200ll ). for example, a 
value of u b ~ 7 x 10 14 Hz is expected for a common X-ray spectrum (of the core and the 
brightest knots, inclu ding; knot A) consistent with a simple power law with a = 1.46 ± 0.05 
([Marshall et al.ll2002t for alternative, see bel ow). For the case of ID = 110 pc (IBiretta et al. 



199ll ) and B m = 4.59 mG (jOwen et al.lll989l ). we see in Figure 2 (top) that the v^- value calls 
for a fluctuation level of b ~ 1 x 10~ 3 , consistent with the assumption of the small- amplitude 
disturbance (§ 2.1); that is why this 6-level has been taken into account for Figure 1 (bottom). 

Let us now evaluate the F u spectral index, which is well defined in the transitive 
frequency range between Uf, and u*. Recalling the canonical index oti[= (p — l)/2] for 
u < Ub, for heuristic purposes, we give the definition of ct2 = — A log i^/ A log z/, where 
A log F v = -axAlogzy + log^ and Alogz/[= log(u*/u b )] = [2fj, - ((3 - 1)/2] log(Z~D/A c ). That 
is, the index is found to be expressed in the form 



a 2 



P 



i + P±l 



(3) 



1 Eq. ((2]) reflects the ratio of the averaged magnetic intensity of small to large filaments; that is, 

1/2 r „ rn _ fl nil -i/n ( p+l )/ 2 . 

~ 4> in a good approximation for A m i n « A c « ID, 



where Amin is the minimum scale of A, which might be compared to the skin depth (§ 3.3) 
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for (i > (f3 — l)/4. Note that equation (j3J) is generic, as it is independent of X C /ID and 
thereby independent of any magnetic field parameters (such as B m and b) . The second term 
of the right-hand side that just represents the steepening of the spectrum is coupled, through 
p, with the first term, to deviate from 1/2 (a 2 = «i + 1/2) for a conventional synchrotron 
model (IKardashevI 1 19621 ; iLongairl Il994l ). As seems more reasonable, the a 2 value can be 
compared to the spectral index measured in a X-ray band. 

In Figure 3, for (3' = 5/3, we plot, for convenience, a 2 ([i = £t S yn) as a function of j3, for 
a given a.\ in the range of 0.5 — 0.8. Clearly, one can find that for ol\ = const, the value 
of a 2 monotonically increases as the value of (3 increases. As denoted in Appendix A, by 
inverting equation (j3J) we can express (3 as a function of a%, a 2 , and (3'. When assigning 
the observational indices of a = 0.67 (in the radio-optical band) and 1.46 (in the Chandra 
regime of 0.5 — 7.5 keV) to a\ and a 2 , respectively, equation (IA1I) with equation flA2[) gives 
(3 = 4.3 for a reliable case of (3' = 5/3 (see also Fig. 3). This /3-value, somewhat larger 
than the value of ~ 2 ex pected in an early stage of the current filamentation instability 



([Montgomery fc Liulll979l ). is amenable to the updated result of a three-dimensional kinetic 
simulat ion that revealed the marked steepening of spectrum in a fully developed nonlinear 
phase (ISilva et al.l 120031 ; see their Fig. 3), reflecting accu mulation of self-orga nized larger 



filaments via the violent coalescence of smaller filaments ( Honda et al.l l2000al ibl) . In fact, 
the observational results indicate that the X-ray spectrum t ends to be steepened when de- 
parting from the active galactic core (Marshall et al.ll2002l ). while the radio-optical index 



likely remains constant (IPerlman et al. 



200 ll ); for the spectral data sets of, e.g., knot HST-1 



[(ai,a 2 ) = (0.71,1.29)], knotD (0.69,1.33), and knot A (0.67,1.57), equation (TAT]) yields 
(3 = 2.4, 2.8, and 6.6, respectively. The gradual increase of (3 (departing from ~ 2 at the in- 
nermost knot) might, more or less, reflect the actual dynamic evo lution of a fi l ament cluster 
involving the coalescence, which has originally been inspected by lOwen et al.l (119891 ). 



The minor component of electrons bound to the small filaments with weak |B| makes a 
negligible contribution to F u in the region of v < v*, but instead, the electrons can escape 
from the small filaments, to be additionally accelerated above 7* (§§ 2.1 and 3.1). The ener- 
getic synchrotron emission from the free (not persistently bound) ele ctrons tends to be diffu- 
sive, and the spectral i ndex, a> 3 , could be close to (3 (independent of p: iToptygin &: Fleishman 
19871 ; lFleishmanll2006l ) . In the case in which around the critical frequency of u* the flux den- 
sity level of the diffuse component is not negligibly smaller than that for the coherent one, 
e.g. for «3 > a 2 the synchrotron spectrum is further steepened above v* to extend the 
gamma-ray tail. More on this topic is discussed in § 3.3 and in Appendices A and B. 



-10- 



3. ENERGY REBOOST TO THE PeV RANGE 
AND GAMMA-RAY EMISSION 

Below we consider the additional DSA of the injected electrons undergoing the nongy- 
rational motion in the filamentary plasma. We evaluate the achievable highest energy (§ 3.1) 
and the upper limit of the photon frequency (§ 3.2), and argue a condition for which the 
high-frequency tail of the gamma-ray synchrotron spectrum could be extended (§ 3.3). 



3.1. Cutoff Lorentz Factor of Energetic- Tail Electrons 



The sufficient condition for which the free electrons can be additionally accelerated in 
situ is that the acceleration time must be shorter than that for the conventional DSA of 
gyrating electrons (for the necessary condition, see §§ 2.1 and 5). In the high-energy regime 
compared to high-/? plasma state, the free electrons spaciously meander, suffering successive 
small deflections by the quasi-static magnetic fields fluctuating in the zeroth order (sustained 
by bulk current filaments). The fields of many fine filaments can be characterized by the 
rms strength of B = ((|B| 2 )) 1//2 , where (• • •) indicates an average on the transverse two- 
dimensional plane of the filaments; an average vector of the local fields likely vanishes in 
the interior of the jet such that | (B) | <C B, whereby the gyroradius can be no longer 
well defined. The key kinetic property in this energetic regime is that the free electrons 
quickly lose coheren ce of the trajectory in the stochastic medium. For the off-resonant 
scattering diffusion (IHonda fc Hondall2005al ). the correct ordering of the mean free path is 
i e ~ e 2 /(e 2 B 2 lD); here note that its anisotropy is small, according to the three-dimensional 
rms deflection. The ratio of the mfp to be expressed as 



4 in _ 3 7 10/iG100pc 



9,e 



10 9 B ID 



It is evident that there is a wide parameter range in w hich 
ordering of B ~ |B| A ^ Ac ~ 10 /iG and ID = 110 pc (IBiretta et al. 



(4) 



For example, for the 



19911 ). the relation of 



< r g,e is safely satisfied in the range of 7 < 10 11 . The shorter mfp (in the shock normal 
direction) is essential for more rapid back-and-forth across the shock, resulting in a shorter 
DSA timescale. 



2 Relating to a majo r difference from the ion kinetics that was examined in the preceding work 
(jHonda fc Hondal l2004al ). one can find the inequality relation of £ e /r g . e <C ^i/ r g,i < £p/ r g,pi where the 
subscripts "£" and "p" indicate arbitrary ions and protons, respectively. In this sense, the quick loss of the 
coherence in electron motion is more pronounced than that in ion motion. 
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If we allow I to take a value less than unity, the modified version of the mean acceleration 
time for the electrons is found to give the scaling of 

^-H^rm'^m 2 - 

The numerical factor of equation ([5]) is correct for f3 = 4.3 and r = 3.2 (expected at knot A), 
although it is merely weakly dependent on these parameters. For, e.g., B ~ 10 pG, ID = 
110 pc, and U = 0.5c, the acceleration timescale (eq. [5]) is, at 7 = 7*(~ 10 9 — 10 10 ), 
favorably shorter than the critical DSA timescale of [r/(r — l)](cr gfi /U 2 ) ~ 1 — 10 yr for 
the Bohm diffusion limit 1^1 Importantly, this ensures that the sufficient condition for further 
energization of the escape electrons with 7 = 7* (and less) is satisfied. In particular, the 
critical energy 7*mc 2 may be associated with the injection energy for the additional DSA 
(§ 2.1), and the conventional DSA of the gyrating electrons in the escape tracks (illustrated 
in Fig. 1, bottom) can be regarded as the "preaccelerator" in the envisagement of the two-step 
DSA. Since the shorter t acc of equation (j3J) can be retained up to the very energetic regime of 
7 ~ 10 11 (in conformity with £ e < r g e ), we expect that two electron populations mediated by 
the injection threshold, namely, the major low-energy and minor high-energy components, 
will be established. Recall here that for the Fermi (type I) fashion, the resulting index p 
(for the power-law dn/de distribution) is independent of the diffusion coefficient, depending 
solely on the scattering center compression ratio r (§ 2.1). Accordingly, for the present model 
supposing ubiquitous Fermi operation in the two energy levels, it is, at least for a high Mach 



number case, exp ected that the power-law distribution likely has a common p ( ]Vainiolll999 



Schlickeiseril2002t ). 



As usual, the additional DSA of the injected electrons operates as long as they are 
confined in the accelerator, where the electron kinetic energy can increase until t acc asymp- 
totically approaches the radiative loss times. This spatiotemporal limit truncates the power- 
law tail of such a dn/de distribution. We note that the spatial limit may be denoted as 
l e ~ D, instead of r 9 e ~ A for the pre-DSA, and as for the synchrotron los s, the dif- 



fuse effects could be nonnegligible (jToptygin &: Fleishman! 119871 ; lFleishmanll2006l ). We here 



take, for convenience, the harmonic average of the synchrotron and IC cooling times to get 



3 Concerning the limitation of i acc that depends on the diffusi on coefficient, there is a general cons traint 
for which the diffusion speed of electrons must be larger than U ( Jokipii 1987 ; Honda fc Honda 2005bl ) . For 
nonrelativistic shocks, indeed, this constraint turns out to be satisfied, as far as the electron mfp in the 
shock's normal direction does not greatly exceed the characteristic length of the density gradient. Also, 
equation ([5|) cannot take unlimitedly smaller values with decreasing 7, because the gyrotrapping effects of 
cold electrons in the local magnetic fields make the approximation of | (B) | <C B worse. The lower limit of 
7 reflects the energy required for escape from the local magnetic traps. 
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ta.cc = £syn[l + ^(^kn/ot)] 1 , where t SJQ is the diffuse synchrotron time derived in Ap- 
pendix B, f is a dimensionless factor (see eq. |B4j ). a = 8nu p h/B 2 , and Qkn /ot is the ratio 



of the Klein- Nishina and Thomson cross sections (e.g., lLongairlll992l . Il994j ). Note that this 
time-balance equation does not explicitly contain B, because both t acc and t syn have a com- 
mon dependence of being proportional to B~ 2 . Solving the spatiotemporal balance equation 
for 7, given ID, U, a, and the characteristic frequency of target photons in the comoving 
frame, v t (involved in o"kn), yields the achievable maximum (cutoff) Lorentz factor of an 
accelerated electron, defined as 7 C = min (7|t acc =t loss , l\i e =i))- 

In Figure 4, for the standard parameter set of ID = 110 pc and U = 0.5c (§ 2), the 
numerical solutions, 7 C , are plotted as a function of v t in the parameter range of 10 1 < a < 
10 6 . In any case, we find that the spatial limit plays a nondominant role in determining 7 C 
(thereby, j c = 7|i acc =t loss ), in contrast to the previous case of determining 7* (§ 2.1). For the 
larger af value, we see that in the Thomson limit of hv t <C mc 2 /7c, where h is the Planck 
constant, the values of 7 C take constants because o"kn/o"t — 1, while for hv t > mc 2 /j c , in 
which okn/ot < 1, the Klein-Nishina effect enhances the value of 7 C . However, 7 C does 
not unlimitedly increase, but eventually saturates, because of the dominance of the diffuse 
synchrotron cooling loss, to take a constant value. For variations of u t , the allowed domain 
of the other variables can be thus written as follows: 

1 > ^ / 100pc\ 1/3 /0.5c\ 2/3 1 fR . 



5 x 10 10 V ID 



Note that the upper bound and af <C 1 reflect the diffuse synchrotron limit and the lower 
bound for af > 1 reflects the Thomson limit. 

In the region in which B < 1 mG, the EC scatterers, including the starlight and dust 
emission (§ 2.1), can lead to a > 1 for 1 < T < 10. As seen in Figure 4, however, in their 
characteristic frequencies of v t ~ 10 14 r Hz (stars) and 3 x 10 13 T Hz (dust) the values of 
7 C turn out to be limited by the diffuse synchrotron cooling rather than the IC cooling, for 
the displayed range of a < 10 6 that would cover the conceivable a-values. This is because 
the Klein-Nishina effects prolong the IC timescale beyond the diffuse synchrotron time. 
Although the CMB provides a lower target frequency iv t — 2 x 10 11 T Hz) that tends to 
fade the Klein-Nishina effects, again it will play a negligible role in determining 7 C , because 
a < 1; the SSC process is also ignorable, as t he source arranges a lower w p h and higher v t 



compared to those of the above EC sources (IStawarz et al.ll2003l ). Consequently, for the 



standard value of (3 = 4.3 and another possible value of (5 = 6.6 (at knot A; § 2.2), we 
obtain 7 C ~ 5 x 10 10 (Fig. 4) and 4 x 10 10 , respectively, exceeding 7*, and then the diffuse 
synchrotron timescales are, for example, for B ~ 10 /iG, estimated as t syn (~ t acc ) ~ 2 yr 
(eq. [5]) and 1 yr, respectively. Even if the /-value deviates from unity (/ < 1), the decrease 
of these 7 C - values is not so significant on account of the dependence of 7 C oc Z 1 / 3 . 
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3.2. Inverse Compton Gamma- Ray Cutoff 

The electrons concomitant with the truncated power-law distribution scatter the target 
photons with frequency v t boosting them to higher frequency regions. Of particular interest 
is the achievable highest frequency of the scattered photons, since a significantly large value 
of 7 C is accomplished in the present scenario. This is the subject of truncation of the gamma- 
ray spectrum. 

In Figure 5 we plot the highest (cutoff) frequency of the IC radiation (in the comoving 
frame), h> ic ^ c , as a function of v t in the parameter range of 10 1 < a < 10 6 . The values of ID 
and U are chosen to be the same as those in Figure 4. In the high-energy and high-frequency 
region of hv t 3> mc 2 /(4'-f c ), we find that v\ c ^ c ~ 7 c mc 2 //i, while for hv t <C mc 2 /(4 / y c ), which 
covers the Thomson limit, z/ icc ~ ^ll^t- Note that in a Klein-Nishina regime satisfying 
f c > t sjn , where 7 C takes a value independent of v t (Fig. 4), the former expression of z/ iCjC 
(oc 7 C ) is independent of v t , as is obviously seen in Figure 5. The scaling of this upper limit 
is written in the simple form 

\ ^ 3 / \ 2/3 

6 X 1 30(_^M (JL) Hz. (7) 



100 pc I \0.5c 

Apparently, this is just the case involved in the EC scenarios due to the stellar and dust 
emissions with v t ~ 10 14 r Hz and 3 x 10 13 T Hz, respectively (§ 3.1, Fig. 5), to safely give 
the ordering of v 1Cfi ~ 10 30 Hz for the reasonable parameter range. Again we point out 
the weak / dependence and the B independence of equation ([7]), and that the CMB-based 
EC and SSC mechanisms will not greatly change the z/ ic c -value. As a consequence, we can 
claim that the z/ icc -value multiplied by an appropriate Doppler factor, which is of order 
unity in the present regime, provides a common cutoff of the IC radiation from the source. 
Concerning the observability of the very energetic photons, the absorption in extragalactic 
background light with subsequent pair production (INikishovl Il962l ; iGould fc Schrederl 119661 ) 
will be significant in spite of the proximity. The relevant issues, including applicability to 
more distant radio galaxies, are discussed in § 5. 

In Figure 6, for ID = 110 pc and U = 0.5c, we present the model vF v spectra of the 
possible IC radiation due to the EC target with v t = 10 14 Hz (stars) and, for convenience, the 
synchrotron radiation for the corresponding case of /3 = 4.3 (solid curve), and another case 
with (3 = 6.6 (dotted curve) as a comparison. The synchrotron spectra h ave been best fitted to 



the observational data of near-infrared to optic al (at the knot A shock; iPerlman et al.ll200ll ) 



and X-ray bands (knot A; iMarshall et al.ll2002l ). As for the IC radiation component, we plot 



the intrinsic (source) spectrum for the possible cases of a = 10 2 (dot-dashed curve) and 
10 _1 (dashed curve), as accommodated by the data recently observed by the HESS telescope 
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( IBeilicke et al.l 120051 ) that are considered to reflect a quiescent epoch of the (unresolved) 
gamma-ray source. Note that the chosen a- values will sufficiently reflect the uncertain value 
of u p h > 10~ 9 ergs cm -3 and rms strength of the fields that trap the major electrons emitting 
the synchrotron photons with the frequency ~ i/j, (at which vF v takes the peak). As seen in 
Figure 6, below the break frequency of IC gamma rays, the v dependence of vF v coincides 
with that for the synchrotron continuum [i.e., vF v oc u^ 3 ~ p ^ 2 = z/ - 33 ], and the spectrum 
takes the local maximum value of vF v ~ 4 x 10 10 Jy Hz for a = 10 _1 at the break frequency 
of ~ 10 25 Hz. Above the break, vF v is significantly suppressed due to the Klein-Nishina 
effects, followed by the asymptotic dependence of vF v oc v~^ p ~ 1 '\nv (~ z/~ L2 for a spectral 
fit), and ends up with the cutoff at v- xcc ~ 6 x 10 30 Hz (eq. [7]). We note that changing the 
/?-value to 6.6 leads to a slight shift of the cutoff to z/i C C ~ 5 x 10 30 Hz (not shown in the 
figure). In sub-GeV to GeV regimes, the curve for a = 10 _1 reaches the sensitivity limit 
of the LAT instrument (Large Area Telescope) on board the Gamma-Ray Large Area Space 
Telescope (GLAST). 

In regard to the synchrotron spectra, the cases of (3 = 4.3 and 6.6 reproduce, via 
equation ([3]), for (a\, fl') = (0.67,5/3) , vF u oc u~ 0A6 and z/ ~°- 57 in the transitive frequency 
range that covers the Chandra band ([Marshall et al.ll2002l ). giving ~ 7 x 10 14 and 1 x 
10 15 Hz, respectively; and in turn, for the self-consistently determined 6-levels (about 1 x 10~ 3 
and 2 x 10 -3 , respectively, for the typical ID and B m \ § 2.2), we numerically find u* ~ 
1 x 10 21 Hz (Fig. 2, middle) and 2 x 10 21 Hz, as well as ~ 1 x 10~ 5 (Fig. 2, bottom) and 
3 x 10 -6 , respectively. As seen in Figure 6, for example, for the case of a = 10 _1 , the vF v 
levels of the synchrotron component are, above the frequencies of 7 x 10 19 and 3 x 10 19 Hz 
(for each of the /3-value cases), overcome by those of the IC component, and at v — v*, they 
are lower than the IC level by a factor of about 4 x 10~ 2 and 1 x 10 -2 , respectively. 



3.3. On the Gamma-Ray Tail of the Synchrotron Spectrum 

Regarding the energetic synchrotron radiation above z/*, we are concerned with the 
following two possible channels of the related electron kinetics in actual circumstances: (1) 
the additionally accelerated electrons mainly retain the meandering motion without getting 
trapped in the filaments, and (2) the electrons (with energy e) are transported to be trapped 
in the larger filaments with the size (A) in the range of A c < 2r ge (e) < A (< D) and are 
cooled down through the gyration with r ff e . For case 1, the diffuse synchrotron spectrum 
exhibiting F v ^ oc can appear in a region of v > v*. At v = u*, the flux density ratio of 
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the diffuse to the coherent synchrotron component is found to scale as 



2 x 10 



-15 



1(T 3 pc 



3.3 



10 



7.3 



10 /iG 
B 



3.3 



for (3 = 4.3 (see eq. |B5j for an arbitrary (3). For the case reflected in Figure 1 (bottom) 
and for the diffuse synchrotron limit (Fig. 4), we have 7 c /7* — 7 and B (~ |B|a=a c ) ~ 
5 /iG, whereupon in equation ([H]), the continuity condition of [F v aI Bv)v=v* c = 1 requires a 
considerably small size (A) ~ 3 x 10~ 8 pc (9 x 10 10 cm). Although the (A)-value is much larger 
than the order of magnitude of the skin de pth (~ 10 6 cm for a plasma density of ~ 10" 1 cm -3 ; 
Owen et al.lll990l ; iDi Matteo et al.l 120031 ) . it is still unclear whether a large number of the 
filaments that are much smaller than ID can actually exist or not. For (F v a/ F v ) v=v * <C 1, 
the critical frequency u* can be regarded as a cutoff accompanied by the exponential-like 
decay of F u , and such a case is shown in Figure 6, as an example. 

For case 2, in which the orbit stochasity of the rms deflection becomes incomplete, the 
diffuse synchrotron spectrum will be modified such that the emission of the fundamental 
mode appears to be preferentially pronounced. For an estimation of the effective frequency 
u c , one can employ the simple formula (oc 7 2 |B|), intending to make a comparison with the 
previously suggested cutoff for coherent synchrotron radiation (§ 1). Combining equation ([6]) 
with the formula yields 



1 > 



100 /iG /100 pc 



1 x 10 24 Hz IBI 



A>A C 



ID 



2/3 



0.5c 



4/3 



> 



1 + af) 2 / 3 ' 



(9) 



In equation ([9]), the upper limit (or ar <C 1), achievable in the knot A environment (§ 3.1), 
gives a scaling of v c oc |B| A>Ac (ID) 2 / 3 U 4/3 , and then, for given values of ID = 110 pc and 
U = 0.5c, the expected range is v c > 10 23 Hz for |B| A>Ac > 10 /zG. Because of the large value 
of 7 n , the value of v c , is found to be strikingly larger than the previous values of ~ 10 15 Hz 
( lKeellll988l : iPerez-Fournon et al.lll988l ; iMeisenheimer et al.lll996l : IStiavelli et al.lll997l ) , which 
could be deduced from a simple DSA mo del assuming the resonant scatte ring diffusion of the 
electrons confined in a uniform jet (e.g., iBiermann fc Strittmatterlll987l for a parallel shock 
case). 



4. Variability of the Broadband Spectrum 

To complement the model predictions, we give a discussion of the variability of the 
broadband spectrum, shown in Figure 6. In the low-frequency bands associated with syn- 
chrotron radiation, it is expected that the light curves reflect the spatial inhomogeneity of 
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the interior. The synchrotron emissions from the escaped electrons with Lorentz factors 
above 7* = [eB m W / \2mc 2 )}{\/lD) {1+13 ^ 2 (i.e., above the dotted line in Fig. 1, bottom) will 
exhibit the high-frequency variation with a typical light-crossing time of order ~ A/c. In 
this regime, the highest (cutoff) frequency of the power spectrum might be compared to the 
plasma frequency. On the other hand, the electrons in the quasi-secular trap tracks (the 
shaded regions in Fig. 1, bottom) cannot freely cross the filaments, in which, instead of such 
a variation with a timescale of A/c, the low- frequency variation due to the shock int ermit- 



tency that might reflect the activity of the galactic central engine (IBiretta et al. 



appear, conforming to the 1/ f — 1/ f 2 flicker-noise spectrum (e.g.. lKataoka et al 



19951] could 
200lh . The 



bound electrons that emit a synchrotron photon with frequency v have a Lorentz factor of 
7 = [47rmcz//(3e J B m )] 1 / 2 (A/[D) (1 - /3)/4 . Equating 7(A,v) to 7*(A), we find the solution A (as 
a function of u), which is for convenience labeled as A& (the high-frequency variation relies 
on the spatioscale of A < A&). The quantity c/A& = /& then characterizes the frequency 
that discriminates between the aforementioned high- and low-frequency regions of the power 
spectrum for the v-b&nd light curve. This key frequency, analogous to the bend-over fre- 
quency of the broken power spectrum for blazar light curves (Appendix C), is found to scale 
as 



fb ^ 2 x 10- 



10 i5 Hz 



0.14 



0.43 



100 pc 



0.71 



Hz 



(10) 



lmGy V ID 

for (3 = 4.3 (see eq. |C1| for the generic expression). Note the very weak v dependence 
of ft,, in accordance with the feature of the blazar light curves. For the parameter set of 
B m = 4.59 mG and ID = 110 pc, equation (02) gives f b ~ 4 x 10~ 8 Hz and 2 x 10" 8 Hz for 
v = 10 15 Hz (optical/UV) and 10 18 Hz (X-ray), respectively. We note that for j3 = 6.6, these 
/^-values shift down to 5 x 10~ 9 Hz and 3 x 10~ 9 Hz, respectively. As a consequence, the 
model predicts the characteristic variation with a 1 yr and 10 yr timescale (for f3 = 4.3 and 
6.6, respectively), longer than a day timescale for typical blazar jets (Appendix C; eq. [Clj ). 

Concerning the MeV/ GeV and sub-TeV/TeV gamma-ray light curves, the high-frequency 
variation (reflecting the filament size) cannot appear, as the bound electrons in the quasi- 
secular trap tracks dominantly contribute to Comptonize the target photons. The variation 
could solely reflect the intermittence of the shock accelerator in the flicker-noise regime. In 
the EC model supposing a stationary u p h value, it then appears that the gamma-ray vari- 
ation must correlate with the variation of the synchrotron spectrum in the low-frequency 
(flicker noise) regime below /&. This feature will provide a constra int to explain the data 
of the High Energy Gamma Ray Astronomy (HEGRA) experiments (lAharonian et al.ll2003l ) 
that seem to reflect the fl aring phase of the source (see Fig . 6) . Even though th e other 



scenarios (IBai fc Ledl200ll for the SSC; iProtheroe et al 



synchrotron-proton blazar [SPB] model; IStawarz et al. 



2003 



and Reimer et al. 



2004 for the 



2006al for the EC involving nuclear 
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activity) could also be candidates for explaining the higher flux level and variation, we can 
claim that the present simple model will guarantee the level shown in Figure 6, which is 
compatible with the lower levels recently confirmed by HESS (IBeilicke et al.ll2005T). as we ll 
as consistent with the updated Whipple 10 m telescope upper limit (ILe Bohec et al.l 12004 ) . 



5. DISCUSSION AND CONCLUSIONS 

In the present conception, the lack of the filaments with A < A c leads to the failure of the 
transition injection followed by the additional DSA. Moreover, even if fine filaments exist, 
the large value of (A) lowers the F ud level, as is deduced from equation (jSj) (or eq. |B5j ). It 
is, therefore, important to confirm the presence of such a small scale in the actual jet. In this 
aspect, finding /& and the high-frequency extension of the power spectrum (for the observed 
light curve; § 4) provide the implication of the presence of the small scale of A < A&, and in 
addition, the truncation of the power spectrum just reflects the actual smallest value of A, 
whose lower limit is expected to be on the order of ~ 10 6 cm (§ 3.3). Another feasible method 
is availing of the property according to which emission from the electrons in the stochastic 
medium could be unpolarized. The threshold frequency above which the polarization de- 
creases is estimated as z/ t h = (3/47r)(|eA|/mc 2 ) 2 (|eB|/mc), where B = Vx A, to give the cor- 
responding coherence length that scales as A ~ 2 x 10~ 3 (z/ t h/10 15 Hz) 1 ^ 2 (10 /iG/|B|) 3 ^ 2 pc. 
Combined with the expression of |B| (as a function of A), the scaling is recast to 

u th \0-i4 / 1 mG x ' ' 




^ m (.TwirJ UrJ pc (11 > 

for (3 = 4.3. In particular, for the knot A parameters of ID = 110 pc and B m = 4.59 mG, 
equation (JTTj) reduces to A ~ 0.3(z/ t h/10 15 Hz) 014 pc, simply having a we ak dependence 



on Uth] then, for example, for the optical range of is t h (jPerlman et al.lll999l ). we have A ~ 
(2 — 3) x 10 _1 pc, smaller than A c ~ 2 pc (or 7 pc for (3 = 6.6). The direct measurement of 
A seems to be difficult even by the modern Hubble Sp ace Telescop e with a spatioresolution 



of 0.1"; viz., 7.8 pc at the distance to M87 (~ 16 Mpc; lTonrylll99ll ). while the A c - value may 
be comparable to the spatioresolution. 

In conclusion, we have given a full argument in favor of regarding the M87 jet/knot A as 
a hard X-ray and gamma-ray emitter, owing to the Pevatron electron accelerator comprising 
the two-step DSA based on the filamentary jet model. In a plausible case in which the 
Alfvenic fluctuations superimposed on the local mean magnetic fields (sustained by the bulk 
current filame nts) are in a Kolmogq rov turbulent state, the observed X-ray spectral index of 



1.46 (or 1.57) (Marshall et al.ll2002l ) has been reproduced for the spectral index of filamen- 
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tary turbulence of (3 = 4.3 (6.6), which is amenable to the results of the kinetic simulations 



that solve a long time evolution of fi 



a shock speed of 0.5c (IBiretta et al 



amentation instability. For the typical para meter set of 
1995) , a field in clination angle of 60° (e.g., lOwen et al. 



19891 : ICapetti et al1ll997l : Irleinz fc Begelma^l Il997h. a comp ression ratio of 3.2 (compati- 



Perlman et a 



ble with the radio-optical index of 0.67; 
compared to the jet diameter of 11 pc (IBiretta et al. 



2001|), a filament correlation scale 



199ll ). and the magnetic intensity of 



an outer scale filament of 4.59 mG (jOwen et al.lll989l ). we have self-consistently determined 
the Alfvenic fluctuation level and found that the synchrotron cutoff most likely appears at 
1 x 10 21 Hz (for (3 = 4.3; 2 x 10 21 Hz for (3 = 6.6), where the flux density is, by a factor of 
1 x 10~ 5 (6 x 10~ 6 ), lower than that inferred from the extrapolation of the radio-optical con- 
tinuum. For the optical/UV and X-ray light curves, the bend-over frequencies of their power 
spectra have been estimated as 4 x 10~ 8 Hz (5 x 10 -9 Hz) and 2 x 10~ 8 Hz (3 x 10~ 9 Hz), 
respectively. The additional DSA incorporated with the transition injection turns out to 
require fine filaments with a transverse size smaller than 2 pc (7 pc) to boost the Lorentz 
factor of a free electron to 5 x 10 10 (4 x 10 10 ), independent of the rms magnetic field strength 
of the filaments. The intrinsic cutoff frequency of the conceivable IC (EC and SSC) radiation 
is worked out at 6 x 10 30 Hz (5 x 10 30 Hz), which just reflects the uppermost electron energy 
of ~ 20 PeV. It is expected that the MeV/GeV and sub-TeV/TeV gamma-ray light curves 
solely reflect the intermittence of the shock accelerator in the flicker-noise regime. 

The present conception for the generation of very energetic emission from the large-scale 
jet will constitute a frontier, complementary to the oth er scenarios that pre dict emission from 
the unresolved galactic core (IProtheroe et al.ll2003l and lReimer et al.ll2004l for the SPB model 
as a misaligned BL Lac object). Concerning potential applicability to more distant sources, 
if the similarity of blazar variability to the present one (eq. |10j ; Appendix C) is confirmed 
in future observations, the filamentary jet model could be useful for inferring the detailed 
structure of blazar jets. It then follows that the generic modeling can be used to reinforce 
an interpretation of the intrinsic vF v spectra that have been less well theoretically studied 
so far. Even for the nearby galaxy M87 [z = 0.00427), the cosmic infrared background 
(CIB) will opaque the gamma rays with energies in the range 10 — 100 TeV, although the 
uncertainty of the extragalactic CIB lev e l makes it difficult to ac c urately determin e the 
optical depth (jCoppi Sz Aharonian 1997 ; Protheroe fc Meyer 2000 ; Aharonian 2001 ) and 
thereby the observed cutoff (the expected range is indicated in Fig. 6). By contrast, the 
PeV photons ought to suffer the absorption dominantly by the CMB, so that their mfp can 
be predicted with very high accuracy, provided that the Lorentz invariance remains even in 
the ultrahigh energy region. In any case, the present model spectrum is readily referable to 
the radio-loud sources within ~ 100 Mpc, for which the extragalactic background (CMB, 
CIB) light is almost transparent for the sub-TeV (TeV) photons dominantly carrying the 
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IC gamma-ray flux, whereas for the more distant sources the observed gamma-ray cutoffs 
could shift to a lower frequency region, down to the lower limit that seems to appear around 
~ 10 GeV fcneiske et al.ll2004h . 



Relating to this perspective, the pair-cascade radiation triggered by the photon-photon 
annihilation potentially contributes to raisin g; the extragalactic gamma-ray background level 
( iSreekumar et all Il998l ; IStrong et all 12004 ) . When assuming a stellar emission level, the 
gamma-ray luminosity apparently depends on the rms magnetic field strength; that is, the 
complexity of the filamentary turbulent state (including a spectral index deviating from that 
of classical hydromagnetic turbulence; eq. |Alj ). But the contribution from all Fanaroff- 
Riley (F R) type I radio galaxies cont aining X-ray jets is expected not to exceed the one from 
blazars fjKneiske fc Mannheim! 120051 ) as long as the rms strength is, on average, not much 
smaller than the equipartition value (like the case sho wn in Fig. 6 for a < I P -1 ), along with 
the arguments invoking the luminosity function (see IStawarz et al.l l2006bl for the details). 
Particularly, for the smaller power index of the filamentary turbulence, the rms strength 
of the fields that trap major electrons constituting the synchrotron vF v peak tends to be 
comparable to the uppermost value of ~ 1 mG for the outer scale filament, wherein the 
smaller contribution of at most 1 % is merely guaranteed. These arguments are not affected 
by the shift of the intrinsic cutoff in the PeV range, as the dominant contribution commonly 
owes to the (sub-)TeV photons constituting the local maximum of vF v . To verify the pre- 
dictions, in addition to the survey of hidden sources by future high-sensitivity detectors, 
comprehensive measure ments of the filamentary turbulence and its spectrum as well would 



be of powerful use (e.g., ICarilli &: Barthellll996l ). bringing on some interesting information 



arguably linked to synchrotron spectra (§ 2.2; Fig. 3). Modern laboratory experiments of 
high-intensity la ser-matter interaction, which can reproduce the similar turbulent state (e.g., 
Wei et al.l 120041 ). might be also available for resolving the complexities in detail. 



A. THE EXPRESSION OF THE SPECTRAL INDEX OF 
FILAMENTARY TURBULENCE 

For instruction, here we write down the expression of the power-law spectral index 
of magnetized filamentary turbulence, j3, as a function of the indices «i, a 2 (> «i), and (5 1 
1). Again, note that «i (for v < z/ fe ) and a 2 (for v h < v < is*) could be typically compared 
to the radio-optical and X-ray indices, respectively, and f3' indicates the power-law index of 
the Alfvenic fluctuation disturbing the zeroth-order magnetic field (permeating through the 
current filaments). If we make use of equation ([I]) for the ordinary case of /i = fi syn ((3, /?'), 
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the inversion of equation ([3]) readily gives 



/3(ati,a 2 ,P') 



x(ax, a 2 ) + x a ((3') 



x(ai,a 2 ) - 3x a ((3') 



x(ati, a 2 ) 



13' -I 



(Al) 



(A2) 



(X 2 — Oil 3 — (3' 

Particularly, for a likely case of (3' = 5/3 (reflecting the Kolmogorov turbulence), we have 

X a = 1/2. 

In conjunction with the synchrotron spectrum, in the very high frequency range of 
v > is*, the F v spectrum could have an index of «3 = (3 (§§ 2.2 and 3.3; eq. |B1| ). As is, it 
may be also instructive to clarify what condition would determine the softening (a^ > a 2 ) 
or hardening (a 3 < a 2 ) of the spectrum above u*. The algebraic manipulations lead to the 
Ansatz that the relation (3 > a 2 {(3 < a 2 ) is obtained for the condition a 2 > a\ (a 2 < a 2 ), 
where the critical index a\ can be expressed as a function of o>\ and f3' in the form 



<4 (<*!>/?') 



x b {ai,(3') + ^x 2 b (a u f3>) - 12 ai x a ((3')x c (f3') 

QXaW) 



(A3) 



where the abbreviations are x b = cti(l + 3x a ) + x c and x c = 1 — x a . In the special case that 
satisfies the equality of a 2 = a 2 , we have the relation o f a 2 = (3 leading to a common index 
of a 3 = a 2 . For example, given a>i = 0.67 (at knot A; iPerlman et al.ll200ll ). equation flA3l) 
for (3' = 5/3 provides a\ — 1.3, whereby t he comparison with t he index of a 2 = 1.46 ± 0.05 
to 1.57 ± 0.10 at the corresponding point ([Marshall et al.l 120021 1 leads to a^> a 2 . 



B. 



THE COOLING TIMESCALE AND FLUX DENSITY LEVEL OF THE 
DIFFUSE SYNCHROTRON RADIATION 



On the theoretical basis by iToptygin &: Fleishman! (119871 ) . we derive the cooling time 
of the diffuse synchrotron radiation of additionally accelerated electrons in the magnetized 
filamentary plasma. The emissivity of a single electron by the diffuse synchrotron radiation 
can be expressed as 



2? {(3 - 1) {(3 2 + 7(3 + 8) e 2 w, 



St 



/3(f3 + 2) z ((3 + 3) c (j V w 



7 2 ^o 



;bi) 



for oj (= 2-kv) ^> 7 2 cg> , where u st = eB/(mc) and ujq = 2nc/ (A) (jFleishmarj 120061 ) . In gen- 
eral, the temporal decay of the electron kinetic energy conforms to —(de/dt) = J du)j(u},e), 
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where e = 7?nc 2 (lLongairlll994l ). For convenience, here we set the lower limit of the us integral 



to 7 2 u;o and introduce the variable transformation of uj/^ujq) = £. Then we have 

„2, . \ P 



doj- 



2 2 

7 <4* 



(3-1 



2 2 

7 



:b2) 



Making use of equation (1B2I) . we obtain the temporal evolution equation of e as follows: 



(B3) 



where <tt is the Thomson cross section, u m = B 2 /8n, and 



f{(3) = (6tt)- 



g(g+ 2) 2 (/? + 3) 
2^ {(3 2 + 7^ + 8) ' 



The solution of equation (1B3|) is found to be e(t) = e (1 + £/£ 



syny 



where e 



(B4) 



7 mc 



= e(t = 0). Here t syn defines the characteristic decay time for the diffuse synchrotron 
radiation, which can be expressed as t Byn = fi syn , where i syn = 3mc/(4crTM m 7o)- Note that 
the form of t n yn appears to be the same as that for the normal synchrotron cooling time 



( lLongairlll994l ). although u m now includes the rms field strength B, instead of a mean field 
strength. For values of, e.g., (3 = 4.3 and 6.6, equation (1B4|) gives f = 3.1 x 10~ 3 and 
1.4 x 10~ 3 , respectively. 

In terms of the diffuse synchrotron radiation spectrum of Jd(u>) = J j(u,e)ne~ p de, we 
pay attention to the ratio to the ordinary synchrotron spectrum [e.g., eq. [18.49] for J(uj) in 
Longairlll994l |. Particularly, at uj* = 2-kv* the ratio can be written in the following form: 



J(u) 



2-p/ 2 (2/3) 



/Ml/2) 



09 



27irg(p) (2/3 + p + 1) f(J3) 



u,d 
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7r 

7* 



2/3-p+l 



(B5) 



where the dimensionless function q (p) including g amma functions is identical with the a(p) of 
equation (18.49) and Table 18.2 in iLongairl (11994J ). It is noted that in the derivation of equa- 
tion ( 1B5I) . the relation of uj st /uj* = (2/3)(7*)~ 2 , which is valid for the ordering B ~ |B|a=a c , 
has been used. For a value of (3 that is not small, the factor {uj st / ujq)^^' 1 ^ can be consid- 
erably small (because typically uj st /uj ^> 1), such that the flux density ratio, (F u ^ d / F u ) u=u *, 
is much smaller than unity (see eq. [8] for (3 = 4.3). By contrast, for a not-so-small value of 
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equation (1B5I) that might be realized at inner knots (such as knot D) with a smaller value of 
(3 ~ 2 — 3 (§ 2.2), the energetic tail of the diffuse synchrotron radiation can appear above 
v*, to exhibit F U)d oc u^ 13 . 



C. THE BEND-OVER FREQUENCY OF THE VARIATION POWER 
SPECTRUM FOR BLAZAR GEOMETRY 



We check the validity of equation (TlOj) by comparing with the corresponding frequency 
of the variation power spectrum for blazar jet geometry with a narrow viewing angle. This is 
associated with (1) application of the filamentary jet model to blazar jets and (2) verificatio n 
of the idea regarding the FR I radio jet as a misaligned BL Lac object (jTsvetanov et al.lll998l ); 
that is, unifying astrophysical jets. If we take the relativistic beaming of the filamentary 
bulk flow into account, the synchrotron frequency is boosted by the beaming factor 5 (i.e., 
v oc #7 2 |B| in the laboratory frame), such that the expression of 7(A, v) denoted in § 4 should 
be also corrected, to be multiplied by <5~ 1//2 . Note that the expression of t£(A) (derived from 
the spatiolimit) is unchanged. For electron acceleration by fully relativistic shocks (with 
the Lorentz factor T), we expect that a strong anisotropy of the momentum distribution 
function is established, such that the accelerated electrons are preferentially directed in the 
narrow cone angle o f order ~ 1/T (in the b ulk rest frame), regardless of the complexity of 
turbulent scatterers (lAchterberg et al.ll200ll ). This implies that the transverse light-crossing 
will be delayed by a factor of T. Therefore, the bend-over frequency in the laboratory 
frame may be redefined as f\ = 5[c/(Y\ b )], where A& denotes the solution A of the equality 
8~ 1 ' 2r y(X, v) = 7* (A) (see also § 4). As a consequence, we find that can be expressed as 



fb 



cS 

Jdt 



-(-) 

Air Kmc/ 



SB r , 



v 



2/(3/3+1) 



eB m lD 
2mc 2 



4/(3/3+1) 



(CI) 



Particularly, setting 8 = Y = 1 and j3 = 4.3 yields the scaling of equation ( flOl) . 



For the bulk Lorentz factor of < T, we have the relation of 5 < T for blazar geometry, 
and then the scaling of equation (1C1I) . e.g., for (3 = 4.3, appears to be equation (fTOj) merely 
multiplied by the factors of 5 /Y (< 1) and 5 2 ^ 3/3+1 ^ = 5 a14 . For the typical blazar parame- 
ters of B m = 10 _1 G, ID = 10~ 2 pc, and 5 = 10, we obtain the value of fb < 5 x 10~ 5 Hz 
for v = 10 18 Hz (X-ray light curve), whic h is reasonably cons istent with the observational 
data (e.g., /[,~3x 1CT 5 Hz for Mrk501; IKataoka et al.ll200ll ). The observational trend of 
the v independence of fb can be also inferred from equation fIClj) for conceivable /3- values. 
According to these, the use of equation (TlOl) will be justified for the concerned case of a not-so- 
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narrow viewing angle and a weak relativistic shock (IBiretta et al.lll995l ; lLobanov et al.ll2003l ). 
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FIGURE CAPTIONS 



Fig. 1.— Lorentz factor of an electron achievable for a conventional DSA 7*, vs. the 
transverse size of a magnetized filament A for a fixed b and some given valu es of B m (as labeled ; 



top) and vice versa (middle), and for b = 10~ 3 and B m = 4.59 mG ( bottom: lOwen et al.lll989l ). 
indicating the regimes where 7* is determined by the spatial (thick dotted line) and temporal 
limits (synchrotron loss; thick solid curve). Here b is the energy density ratio of the Alfvenic 
fluctuation/local mean magnetic field, B m is the magnetic intensity of the filament with A 
reaching a c orrelation length com pared to the jet diameter, D = 110 pc at knot A (thin long- 



dashed line; iBiretta et al.lll99ll ). and more on the DSA parameters (at the knot A region) 
is given in the text. The axes are common in the figures. Note that the solid curve in the 
top panel is the same as that in the middle panel. In the bottom panel, the thin dot-dashed 
curve indicates the (pseudo) limit due to the Thomson scattering, which is removed by the 
Klein-Nishina effect, and the accelerator is active in the shaded domain, wherein the solid 
arrows represent the pre-DSA of the electrons bound to the filaments with A < A c , leading to 
the escape (bound-free transition) entailed in the injection of additional DSA, and the dotted 
arrows represent the energization that is knocked down by the synchrotron loss before the 
escape, namely, the tracks of secular traps in the filaments with A > A c . [See the electronic 
edition of the Journal for a color version of this figure.] 

Fig. 2. — Dependencies of the break (top), the critical frequency v* (middle), and the 
flux density reduction factor (defined in eq. [2]; bottom) on B m , for b = 10~ 4 — 10~ 2 (as 
labeled) and the correlation length of filaments compared to D (= 110 pc; thick curves) and 
0.1D (thin curves). Here the other DSA parameters are the same as those for Fig. 1. The 
horizontal axis is common in the figures. Note that in the middle panel, the thick solid and 
dot-dashed lines overlap the thin ones, and in the bottom panel, the thick dot-dashed line 
overlaps the thin dotted line. 

Fig. 3. — Transitive synchrotron spectral index a 2 (eq- [3]) vs. the power-law spectral index 
of filamentary turbulence j3 for the given canonical index «i as a parameter (as labeled). In 
connection with observed synchrotron spectra, oti and 02 can be typically compared to the 
radio-optical and X-ray spectral indices, respectively. For further explanation, see the text. 

Fig. 4. — Cutoff Lorentz factor of an additionally accelerated electron 7 C vs. the target 
photon frequency v t , for the energy density ratio of target photons/magnetic fields in the 
range a = 10 1 — 10 6 (as labeled). Here the parameters relevant to shock and filamentary 
turbulence are the same as those for Fig. 1 (see also the text). The dotted line indicates 
e = ^Jivti '(me 2 ) = 1: the region of e 1 reflects the Thomson limit, in which the values of 
7 C take constants, while e > 1 indicates the Klein-Nishina regime. An asymptotic value for 
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ar <C 1 (r ~ 3 x 10 3 ) provides the upper bound of 7 C (independent of the rms magnetic 
field strength B), which is referred to as the diffuse synchrotron limit. 

Fig. 5. — Cutoff frequency of the inverse Compton gamma-ray radiation Ui C)C vs. the target 
photon frequency v t for a = 10 1 — 10 6 (as labeled), corresponding to the cases shown in Fig. 4. 
The ranges of characteristic frequency for the candidate external Compton (EC) sources 
of the starlight ("Star"), dust emission ("Dust"), and the cosmic microwa ye background 



("CMB") are also indicated for the bulk Lorentz factor of 1 — 10 (arrows: Stawarz et al. 



20031 ). 



Fig. 6. — Model vF v spectrum of synchrotron radiation for (3 = 4.3 (solid curve) and 6.6 
(dotted curve), and for the former, the source spectrum of IC radiation for the possible 
cases of a = 10" 1 (dashed curve) and 10 -2 (dot-dashed curve) in the knot A environment. 
Here the EC target photon frequency has been supposed to be v t = 10 14 Hz, and the other 
conditions are the same as those for Figs. 1 (bottom) and 4. T he open squares and dia- 
mond i ndicate the data mea sured in the near- infrared to optical (IPerlman et al.ll200ll ) and 
X-ray ([Marshall et al.ll2002l ) bands, respec tively, and the open a nd filled tri angles indicate 



the g amma-ray levels detected by HEGRA (lAharonian et al. 



20051 ). respectively. In addit ion, the EGRET (IReimer et al. 



2003) and HESS ( Beilicke et al. 



20031) and Whipple 10 m tele- 



scope (ILe Bohec et al.l 120041 ) upper limits are indicated. For convenience , the sensitivity 



jerating Cerenkov te lescope (HESS; iHofmannl 120011 ) and future de 



McEnery et al.ll2004l ; New X-Ray Telescope Hard X-Ray Imager and 



limits of a currently o 
tectors (GLAST LAT, 

Soft Gamma-Ray Detector [NeXT HXI and SGD], iTakahashi et al.l 120041 ) are also depicted. 
The intrinsi c cutoff and expected range of the observed cutoff a s a result of the a bsorption 
in the CIB fjCoppi fc Aharonianl 119971 ; iProtheroe &: Meyerl |2000| ; lAharonian! l200ll ) are indi- 
cated by the thick arrow and the thick bar with thin dotted lines, respectively. Note that 
the diffuse synchrotron component that appears in the high-energy gamma-ray region is of 
a considerably lower level. [See the electronic edition of the Journal for a color version of 
this figure.] 
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